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(*_») ■ Parton distributions for the pion are studied in a chiral quark model character- 

^J ' ized by a quark propagator for which a spectral representation is assumed. Elec- 

Cn ' tromagnetic and chiral symmetry constraints are imposed through the relevant 

Ward-Takahashi identities for flavoured vertex functions. Finiteness of the theory, 
'~{ , requires the spectral function to be non-positive definite. Straightforward calcu- 

lation yields the result that the pion structure function becomes one in the chiral 
limit, regardless of the details of the spectral function. LO and NLO evolution 
provide a satisfactory description of phcnomenological parametcrizations of the 
valence distribution functions but fails to describe gluon and sea distributions. 
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1 Introduction 
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r — ' Deep inelastic scattering (DIS) is regarded as one of the traditional ways to 

unveil the quark substructure of hadrons El. The hadronic tensor for electro- 
production on a (unpolarized) hadron is defined as 

^ ; w„(p, q ) = i- f^e^ip\[jjr{ajr(o)]\p) a) 

Ph, j 

Q_i| where p is the hadron momentum, J^ m (£) = q(£,)Ql^l(0 the electromagnetic 

current with Q = diag(e u , e^, e s ) the quark charge matrix and q the photon 
momentum transfer. Introducing the relativistically invariant Euclidean mo- 
mentum transfer Q 2 = —q 2 and the Bjorken variable x = Q 2 /2p • q, in the 
K> , Bjorken limit ( Q 2 — > oo with x fixed), Eq.(Q) can be expressed in terms of 

5^ ' the quark and antiquark distribution functions, qi(x) and q~%(x) 



W^F(x) 



^^ — -f^-YxJK^-Tx 



(2) 



with F(x) = ^ li ef [qi(x) + q~i(x)} /2 as a consequence of Lorentz gauge in- 
variance, scaling and the underlying spin 1/2 nature of quarks. Logarithmic 
scaling violations due to perturbative QCD radiative corrections 13, relate in 
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a linear fashion the structure functions at a given reference scale, Qg, to the 
scale of interest, Q 2 , 

F(x,Q 2 ) = U(Q 2 ,Q 2 )F(x,Q 2 ). (3) 

Here, U(Q 2 ,Qq) is a linear matrix operator, fulfilling renormalization group 
properties U(Ql,Q$)U(Q%,Ql) = U{Q\,Q\) and U{Q 2 ,Q 2 ) = 1 and which 
can be evaluated for high Q 2 . Thus, one is left with the calculation of structure 
functions at a certain scale Qo (the so called initial condition, F(x, Qq)) and 
the assumption that the scale is high enough to make perturbative matching 
meaningful. 

In spite of the wealth of experimental and theoretical information BH'Ufl'Q, 
first principle calculations of the initial condition, F(x,Q 2 ) ), remain, so far, 
elusive. Thus, theoretical quark models of hadrons and quark-hadron dual- 
ity assumptions are used. This obviously raises embarrassing questions on 
confinement, because the intermediate states in the commutator of Eq.dfi) 
should be physical hadron states and, instead, one uses the virtual Compton 
amplitude □ 

T^{p,q) =%Jd i Ze i «*(p\T{JZ a (£)J? a (0)}\p) (4) 

with W^v = l/(27r)ImT M i,. Introducing light-cone (LC) coordinates x = 
{x + ,x~ ,x±) with x ± — x° ± x 3 one gets the following expression for the 
initial condition 

' f ' Ur ~' (21 - " A Vx<P,*) (5) 



F{x,Qi) = -— / " tr 

where the forward quark-target scattering amplitude is defined 

X (p,k) = -i fd^e^ k (p\T{q(Oq(0)}\p) (6) 



x(Pi k) corresponds to the unamputated vertex. Eq.(@) holds under the as- 
sumption of scaling and finiteness in the Bjorken limit. 

Among the hadrons, pseudoscalar mesons and, in particular the pion, ap- 
pear as the simplest and best understood states from a theoretical viewpoint, 
since most of its properties can be explained in terms of spontaneous chiral 
symmetry breaking. Its bound state properties are not expected to depend 
crucially on confinement, at least for small excitation energies. Furthermore, 
the pion_s±ructurc functions have been deduced from phenomenological QCD 
studies W and there also exist lattice calculations of the first few moments 
E3. Thus, we do not expect to achieve a better theoretical understanding of 
a DIS process than in the case of the pion. Hence, structure functions of 
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.the |P' l |V n |J 1 P| vl M hfTV [Computed several times in different quark-loop models 
EiHiS'E^lEJiia'E^Oiil'tj. The issue of implementing gauge invariance and hence 
proper normalization of F(x, Qq) is problematic and therefore requires some 
special care. 

2 Quark propagator and Ward-Takahashi identities 

I present here another model cB for pion structure based in a spectral represen- 
tation of the quark propagator and chiral Ward-Takahashi identities (WTI) 
for the flavoured vertex functions, deduced from conservation of the vector 
current (CVC) and partial conservation of the axial current (PCAC) 

J& a (x) = q(xh^q(x) «'» = (7) 

J% a (x) = q(x)^j 5 y?(x) d^,J^' a (x) = M q(x)iT a j 5 q(x) (8) 

with Mq the average up and down current quark masses. 

We assume a spectral representation for the quark propagator 

S(p)= [dw-?^^-* [d 4 xe lpx (0\T{q{x)q(0)}\0) (9) 

J p-w J 

The requirement of a physical Hilbert space would imply the positivity of 
the spectral function p(w), but we are not assuming this here. Predictions 
based on Eq. (H) require a specific form of p(w) as input, but it turns out 123 
that in the chiral limit, Mg — > in Eq. m), many properties of interest are 
independent of p(w). Thus we will focus on those observables and consider 
the chiral limit. The vector and axial unamputated vertex functions satisfy 
the WTI 

( P / -p),S(p')T^(p',p)S(p) = ^S( P ) 5(p')y (10) 

(p> p)»S(p')r% a (j>',p)S(j>) = 5(p')y75 + 75 y S(p) (11) 

The gauge technique, in the form introduced in QED c3 consists of writing a 
solution for the vector and axial vertices in the form 

1 nTa 1 

7 M 

p 1 ' — w 2 p — w 

p 1 ' — w \ q 2 J 2 j> — w 

with q = p' — p. An interesting consequence of the axial Ward identity ( J13| ) 
is the presence of a massless pseudoscalar pole, which is identified with the 
pion. 



S{p')T^{p',p)S{p) = / d wp ( w )-j_ -y^—L- (12) 

S{p')T^{p',p)S{p) = / dwp(w)-^- ( r " ^ ) 75^— (1-' J .) 
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3 Properties of the spectral function 

Some properties on the spectral function p(w) can be derived by studying 
specific processes and comparing with known QCD results. For instance, 
vacuum polarization can be computed closing the vector-quark-quark vertex 
in Eq. © 

t J d^xe-^iOlT {J^(x)jr(0)} |0) = (-g MV + ^\ II om (p 2 ) (14) 

The e + e _ — ■+ hadrons cross-section is proportional to the imaginary part of 
the vacuum charge polarization operator. Asymptotically one has 

a(e + e~ — ► hadrons) — > — — zZ e ? / dwp(w) (15) 

Thus, the proper QCD asymptotic result is obtained if 

dwp(w) = 1 (16) 



This condition is equivalent to impose that lim p ^ 00 ^S'(p) = 1. 

The weak pion decay constant can be computed from the axial-axial cor- 
relation function, since 

(0 \J^(x)\ n(p)) = i/vP^e** (17) 

and then, inserting a complete set of eigenstates between the currents, 

d 4 xe-^ x (0\T {j% a {x)r A \0)} |0) = ifj a X4- + ■■■ (18) 

where the dots indicate regular pieces in the limit p 2 — > 0. Again, closing the 
axial-quark-quark vertex in Eq. ( J13| ) to make a quark loop, and going to the 
pion pole, p 2 — -> 0, one gets 

L = WfJ dW ^-iJ (27r) 4 (p2 _ W 2 ? (19) 

The momentum integral is logarithmically divergent. Thus, if 

dww 2 p(w) = then f 2 = —-, ^ / dww 2 (- log w 2 )p(w) (20) 

(4tt)^ J 

As we see, p(w) cannot be a positive definite function, regardless of its support. 
f 2 is positive if p(w) is positive around w = 0. Notice also that the first 
condition of Eq. ( |20|) implies that the value of /„. does not depend on any 
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particular scale used to make the logarithm dimensionless. If p(w) had only 
one zero, the scale of /„■ would be set by it. 

A direct application of the previous conditions (ttq) and (EG) can be found 
in the calculation of the pion electromagnetic form factor F° m (q 2 ) which com- 
plete evaluation requires an explicitly knowledge of the spectral function p(w). 
Besides a proper normalization, F° m (0) = 1, the mean squared radius becomes 



ldF 

6 dq 



(rl) = ~^s . - = T^H I dwp(w) (21) 



N c 
q z=o 4tt 2 /^ 

which is the value obtained using the unregularized quark loop if the nor- 
malization condition ( |16|) is used. Thus, the pion has a finite size in this 
framework. 

The amplitude for the (anomalous) neutral pion decay 7r°(p) — ► 7(pi, p) + 
7(j>2, v) can also be computed from the triangle diagram yielding for on-shell 
massless pions p\ — p\ = 0, the result 

n r 

M fJ/U (p 1 ,p 2 ) = % £^ a /3P2Pi / dwp(w) (22) 

which, again, coincides with the standard one if Eq.(nq) is used. 

4 Pion Structure Function 

An advantage of our method is that calculations can be directly done in 
Minkowski space and hence we may directly profit from LC coordinates. For 
the evaluation of the structure function of the pion the relevant vertex is 
constructed from two axial currents and is defined 

(27r) i S(p' + q '-p-q) rf$ v *(jf,Ap,q) (23) 

,r,^x'd 4 y'd 4 y (0\T{j^ a (x)jf(x')q(y)q(ij')}\0)e^'- x ' + P'-y'-^-^ 
and fulfills the WTI 
*rt'rV, fl'jP. l) = e bac S{p')T v v c (p\ p )S(p) (24) 

+ y75^(p' ~ q)n\p' - q,p)S(p) + S 
Up to transverse pieces one gets the solution 

h f 'If WQ Q^ 

XaT' (P .« ;P.«) = / dwp(w)-jj <^ /2 o $ab (25) 

J p — w I q' q z 



~ lb S(p' - q)Y${jt - q,p)S(p) + S(p')T" A b (p',p + q)S(p + q)^ 7 5 



2wq l2 \ n 1 / 2wq^ 



> + — 72- 75 T TT- r-— 7m o 



g ' 2 / ' J 2^ + ^- W V 9 2 / 2 J^'-w 



75 y + crossed}— ■ 
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At the pion poles, q , q' — ► 0, we get the irq — > irq unamputated amplitude 
in the forward direction (a = q') 



u nf i 



X%7'"(P, i;P, l) - q -jrflx ba {p, q) (26) 

and direct use of Eq. flq) together with the conditions ( jig ) and ( J2Cj| ) yields the 
initial condition for the quark distributions in 7r + , 

U7r (x, q§) = 4(1 _ X} qd = e( x )e{\ - x) (27) 

This result has been derived in the chiral limit, is independent of the spectral 
function p(w) and has proper support and normalization. By construction, 
Eq. ( p7J ) is consistent with chiral symmetry, although is not necessarily a conse- 
quence of it. The result has previously been obtained by several means within 
the NJL model either using Pauli-Villars regularization l 13 H 17 1 on the virtual 
Compton amplitude (0) or imposing a transverse cut-off £3 upon the quark- 
target amplitude (||). Within the present approach it can be shown that the 
forward Compton scatterips; amplitude (0) in the Bjorken limit yields the fi- 
nite result presented here Ej based on EqT^). To really appreciate this point, 
let us mention that Eq.(p7|) disagrees with other NJL calculations, due to 
the use-jaf different regularizations. If Eq.(^) is used with a four-dimensional 
cut-off liil or Eq.(|[) is used with Lepage-Brodsky regularizaticp.llj, different 

a apes for the quark distributions are obtained. The null-plane L3, NJL model 
and spectator model 113 calculations also produce different results. In all 
cases, the use of momentum dependent form factors or non-gauge invariant 
regularizations make the connection between Eq. <Q) and Eq. (JsJ) doubtful and, 
furthermore, spoil normalization. The rosnlts based on a quark loop with 
momentum dependent quark masses l 18 H 19 l seem to produce a non-constant 
distribution. 



5 QCD evolution 

Eq. (pTJ) can also be understood t3 in terms of phase space arguments and point 
couplings (i.e., constant matrix elements) which for N constituents gives 

xV(x, Ql) = N{N - l)x(l - x) N - 2 (xV) = / dxV (x, Ql)x = 1 (28) 

Of course, it is tempting to look at the nucleon case N = 3, although we ob- 
viously expect mass corrections to be more important both in the momentum 
dependence of the Nqqq matrix element and in the phase space, and, in addi- 
tion we do not have as much theoretical support as in the pion case. As usual, 
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Figure 1. Valence, gluon and sea distributions in the pion (left) and the nucleon (right) at 
Q 4 = 4GeV compared with phenomenological analysis for the pion SMRS92 7 , GRV99 8 
and for the nucleon GRV95 3 . We take (xV) n = 0.47 and (xV) N = 0.40 at Q 2 = 4GeV 2 . 
The sea and gluon distributions fulfill S(x, Qq) = G(x, Qq) = 0. 



we evolve Eq. (28) at LO and NLO exactly implementing the group properties 
mentioned after Eq.(p3) and employed to look for low energy unitarity limits 
£3. Assuming vanishing initial sea and gluon distributions, the result for the 
dynamically generated valence, sea and gluon distributions for a hadron with 
N = 2 apjl N — 3 constituents compared with known parameterizations for 
the pion Eta and the nucleon u respectively can be seen in Fig p) . The agree- 
ment with the phenomenological valence distributions is good. Actually, since 
normalization of the valence distributions is preserved by evolution and Q\ is 
determined by constraining the valence momentum fraction at Q 2 — 4GeV 
to the phenomenological ones, (xV)^ — 0.47 and (xV)n = 0.40 the bulk of 
the distribution is reproduced. Sea and gluon distributions do not obey such 
strong constraints and provide more information regarding the quality of a 
model. As we see, they are not as good. Nevertheless, a more definite state- 
ment might be made ifjaiicertaintics in the distributions were considered on 
the QCD evolution side O. The estimate of systematic errors in the model, i.e. 
corrections to the initial condition, remains yet a difficult and open problem. 
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